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To what extent is the entropy-area law universal ? 

Multi-horizon and multi-temperature spacetime may break the 

entropy- area law 
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It seems to be a common understanding at present that, once event horizons are in 
thermal equihbrium, the entropy-area law holds inevitably. However no rigorous verification 
is given to such a very strong universality of the law in multi-horizon spacetimes. Then, 

O based on thermodynamically consistent and rigorous discussion, this paper suggests an evi- 
dence of breakdown of entropy-area law for horizons in Schwarzschild-de Sitter spacetime, in 
which the temperatures of the horizons are different. The outline is as follows: We construct 
carefully two thermal equilibrium systems individually for black hole event horizon (BEH) 
and cosmological event horizon (CEH), for which the Euclidean action method is applicable. 
The integration constant (subtraction term) in Euclidean action is determined with referring 
O , to Schwarzschild and de Sitter canonical ensembles. The free energies of the two thermal 

systems are functions of three independent state variables, and we find a similarity of our 
5^ ' two thermal systems with the magnetized gas in laboratory, which gives us a physical under- 

{^X), standing of the necessity of three independent state variables. Then, via the thermodynamic 

consistency with three independent state variables, the breakdown of entropy-area law for 
CEH is suggested. The validity of the law for BEH can not be judged, but we clarify the key 
issue for BEH's entropy. Finally we make comments which may suggest the breakdown of 
entropy-area law for BEH, and also propose two discussions; one of them is on the quantum 
statistics of underlying quantum gravity, and another is on the SdS black hole evaporation 
from the point of view of non-equilibrium thermodynamics. 

CN ' Subject Index: 451, 454 



§1. Introduction 



Entropy-area law, which claims the equilibrium entropy of event horizon is equal 
to one-quarter of its spatial area in Planck units, ^^"^^ is the equation of state of 
^ I the event horizon in thermal equilibrium. This law has already been verified for 

' spacetimes possessing a single event horizon. Then we may naively expect that 

the entropy- area law holds also for multi- horizon spacetimes, once every horizon is 
individually in thermal equilibrium. This expectation is equivalent to consider that 
the thermal equilibrium of each horizon is the necessary and sufficient condition to 
ensure the entropy-area law for each horizon. However this expectation has not been 
rigorously verified in multi-horizon spacetimes. (Comments on existing researches on 
Schwarzschild-de Sitter spacetime will be given in fourth, fifth and sixth paragraphs 
in this section. Please wait for those paragraphs if the reader cares about existing 
researches on multi-horizon spacetimes.) At present, there remains the possibility 
that the thermal equilibrium may be simply the necessary condition of entropy-area 
law. If we find an example that some event horizon does not satisfy the entropy-area 
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law even when it is in thermal equilibrium, then we recognize the thermal equilibrium 
as simply the necessary condition of the entropy-area law. 

We can consider a situation in which the entropy-area law may break down in 
multi- horizon spacetime *\ To explain it, it is necessary to distinguish thermody- 
namic state of each horizon and that of the total system (multi-horizon spacetime) 
composed of several horizons. Even when every horizon in a multi-horizon spacetime 
is in an equilibrium state individually, the total system composed of several horizons 
is never in any equilibrium state if the equilibrium state of one horizon is different 
from that of the other horizon. For example, if the temperatures of horizons in a 
multi-horizon spacetime are different from each other, then a net energy flow arises 
from a high temperature horizon to a low temperature one. Such multi-horizon space- 
time can not be understood to be in any equilibrium state, since, exactly speaking, 
no energy flow arises in thermal equilibrium states. The thermal equilibrium of total 
system (multi-horizon spacetime) is realized if and only if the temperatures of all 
composite horizons arc equal. Therefore, when the temperatures of horizons arc not 
equal, the multi- horizon spacetime should be understood as it is in a non- equilibrium 
state. Here it should be noticed that, generally in non-equilibrium physics, once the 
system under consideration comes in a non-equilibrium state, the equation of state 
for non-equilibrium case takes different form in comparison with that for equilibrium 
case. Especially the non-equilibrium entropy deviates from the equilibrium entropy 
(when a non-equilibrium entropy is well defined). Indeed, although a quite gen- 
eral formulation of non-equilibrium thermodynamics remains unknown at present, 
the difference of non-equilibrium entropy from equilibrium one is already revealed 
for some restricted class of non-equilibrium systems.®-*' Hence, for multi-horizon 
spacetimes composed of horizons of different temperatures, it seems to be reasonable 
to expect the breakdown of cntropy-arca law. However, since a "non-equilibrium 
thermodynamics" applicable to multi-horizon spacetime has not been constructed 
at present, we need to make use of "equilibrium thermodynamics" to investigate 
thermodynamic properties of multi-horizon spacetimes. 

Motivated by the above consideration, this paper treats Schwarzschild-de Sitter 
(SdS) spacetime as the representative of multi-horizon spacetimes. We construct two 
thermal equilibrium systems in SdS spacetime; one of them is for black hole event 
horizon (BEH) and another is for cosmological event horizon (CEH). Note that, since 
the temperature of BEH is always higher than that of CEH in SdS spacetime, we 
need a good way to obtain thermal equilibrium systems of BEH and CEH. As will be 
explained in detail in Sec.2, we will adopt the same way of constructing two thermal 
equilibrium systems as Gibbons and Hawking have used in calculating the Hawking 
temperatures of BEH and CEH;^*^) it is to introduce a thin wall between BEH and 
CEH which reflects perfectly the Hawking radiation coming from the horizons. The 
region enclosed by the wall and BEH (CEH) settles down to a thermal equilibrium 
state, and we obtain two thermal equilibrium systems separated by the perfectly 



*' AH discussions in this paper are based on the ordinary general relativity. The other modified 
theories of gravity are not considered. Even if there is a breakdown of entropy-area law due to exotic 
fields of modified theory, such a breakdown in modified theory is out of the scope of this paper. 
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reflecting wall. Then we will examine the entropy-area law for the two thermal 
equilibrium systems individually. (Although we are motivated by a non-equilibrium 
thermodynamic consideration in previous paragraph, the whole analysis in this pa- 
per is based on equilibrium thermodynamics and we discuss the two equilibrium 
thermodynamics for BEH and CEH individually.) As will be explained in Sec. 2, our 
two thermal systems are treated in the canonical ensemble to obtain the free energies 
of BEH and CEH. Hence we will make use of the Euclidean action method which 
is regarded as one technique to obtain the partition function of canonical ensemble 
of quantum gravity.^^^ (See Appendix A of this paper or of the previous paper. 
Then we will find that the free energies are functions of three independent state 
variables. The existence of three independent state variables for SdS spacctimc was 
not recognized in existing works on multi-horizon spacctimcs.-'^'^)^ -^^^^^^^ But in this 
paper, thermodynamically rigorous analysis with three independent state variables 
will suggest a reasonable evidence of breakdown of entropy-area law for CEH. The 
validity of the law for BEH will not be judged, but wc will clarify the key issue for 
BEH's entropy. These results imply that the thermal equilibrium of each horizon 
may not be the necessary and sufficient condition but simply the necessary condi- 
tion of entropy-area law. The necessary and sufficient condition of the law may be 
implied via some existing works as noted below: 

Let us note that some proposals for thermodynamics of BEH and CEH in SdS 
spacetime are already given to a case with some special matter fields and for an 
extreme case with magnetic/electric charge. These examples are artificial to vanish 
the temperature difference of BEH and CEH, and show that the entropy-area law 
holds for SdS spacetime if the temperatures of BEH and CEH are equal. However, 
in this paper, we consider a more general case which is not extremal and does not 
include artificial matter fields. In all analyses in this paper, the temperatures of 
BEH and CEH remain different and the discussions in those examples^^^ can not be 
applied. If we find the breakdown of entropy-area law for the case that horizons have 
different temperatures, then it is suggested that the necessary and sufficient condition 
of entropy-area law is the thermal equilibrium of the total system composed of several 
horizons in which the net energy flow among horizons disappears. 

Next let us make comments on the case that horizons have different tempera- 
tures. The construction of SdS thermodynamics with leaving horizon temperatures 
different has already been tried in some existing works. Those works as- 
sume some geometrical conserved quantities to be state variables of SdS spacetime, 
and derive the so-called mass formula which is simply a geometrical relation and 
looks similar to the first law of black hole thermodynamics. However, the thermo- 
dynamic consistency has not been confirmed in those works. Here "thermodynamic 
consistency" means that the state variables satisfy not only the four laws of thermo- 
dynamics but also the appropriate differential relations; for example, the differential 
of free energy Fo with respect to temperature To is equivalent to the minus of entropy, 
So = —dFo/dTo *^ . If one asserts some theoretical framework to be a "thermody- 



*' There are many other similar differential relations in thermodynamics. Those relations are 
the ones required in the "thermodynamic consistency", and necessary to understand thermody- 
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namics", the framework must satisfy the thermodynamic consistency. Therefore, 
exactly speaking, it remains unclear whether those existing works^*^^' ""^^^"-^^^ are ap- 
propriate as "thermodynamics" *•*. On the other hand, it seems to be preferable 
that the number of assumptions for thermodynamic formulations of BEH and CEH 
is as small as possible. In order to introduce the minimal set of assumptions which 
preserves thermodynamic consistency, we will refer to Schwarzschild thermodynam- 
ics formulated by York^^ (and also refer to de Sitter thermodynamics^^ which is also 
based on York^^). 

Furthermore, motivated by the dS/CFT correspondence conjecture, some ex- 
isting works^^) focus their attention on the future and past null infinities in SdS 
spacetime (see Fig.l shown in Sec. 2, is the null infinities). Those infinities may 
be appropriate to discuss some geometrical quantities. However, as implied by the 
causal structure of SdS spacetime, the future null infinity seems to be inappropri- 
ate to discuss thermodynamic properties of BEH and CEH, because any observer 
near future null infinity (not near the future temporal infinity i'^) can not "access" 
BEH **). 

Hence, contrary to the existing works, the analysis in this paper is based on the 
following two points: 

• As will be explained precisely in next section, we focus our attention on the 
region enclosed by BEH and CEH (not on null infinity) in SdS spacetime as the 
object of thermodynamic interests. 

• We have a high regard to the "thermodynamic consistency" preserved by the 
minimal set of assumptions without referring to some geometrical conserved 
quantities and dS/CFT correspondence. 

Then, as the result of these two points, a suggestion of breakdown of entropy-area 
law will be obtained. 

Here let us emphasize that, in next section, we exhibit explicitly the assump- 
tions on which our discussion is based. We think readers can judge the approval or 
disapproval to every part of our assumptions and analyses. Therefore, even if some 
part of our discussion and analysis is not acceptable for some reader, we hope this 
paper can propose one possible issue about the universality of entropy-area law. 

This paper is organized as follows: Sec. 2 introduces the minimal set of assump- 
tions in which, with referring to Schwarzschild and de Sitter canonical ensembles,^)' 
the two thermal equilibrium systems for BEH and CEH are constructed, and the spe- 
cial role of the cosmological constant is also pointed out. Sec.3 is for the calculation 
of Euclidean actions of the two thermal equilibrium systems. Sec.4 discusses thermo- 
dynamics for BEH and examines the entropy-area law for BEH. In that section, the 



namic properties of the system under consideration; e.g. phase transition, thermal and mechanical 
stabilities, and equations of states. 

*-* Those existing works^"-*' i^^^^^) preserve/assume the entropy-arca law without confirming ther- 
modynamic consistency. Hence, if the breakdown of entropy-area law is concluded via the thermo- 
dynamic consistency, those existing works can not be regarded as "thermodynamic" theory. 

**^ The observer going towards the future temporal infinity in Fig.l can "access" BEH, since 
the BEH becomes a boundary of the causally connected region of that observer (the region I in 
Fig.l). 
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validity of the law for BEH can not be judged, but the key issue for BEH's entropy 
is clarified. Sec.5 proposes a reasonable evidence of the breakdown of entropy-area 
law for CEH. Finally Sec. 6 is for summary and discussions, in which we make some 
physical comments which may suggest the breakdown of entropy-area law for BEH 
without rigorous verification, and propose two discussions; one of them is on the 
quantum statistics of underlying quantum gravity, and another is on the SdS black 
hole evaporation process from the point of view of non-equilibrium thermodynam- 
ics. And four appendices support the main text of this paper: Appendix A and B 
summarize, respectively, the Euclidean action method and the essence of York's 
Schwarzschild thermodynamic, Appendix C exhibits useful differential formulas for 
calculations of thermodynamic state variables, and Appendix D analyzes the Nariai 
limit (extremal limit) of our SdS thermodynamics. Throughout this paper we use 
the Planck units, c = G = h = kB = ^- 

§2. Minimal set of assumptions 

2.1. Preliminary 

Let us dare to start this section with the discussion given already in de Sitter 
canonical ensemble,^) because this discussion is conceptually essential for thermody- 
namic consistency. 

The aim of this section becomes clear by considering the relation between ther- 
modynamics and statistical mechanics. ^^^'^^^ In statistical mechanics, the partition 
function can not be expressed as a "function of state variables" unless the appropri- 
ate state variables, on which the partition function depends, are specified a priori. 
To understand this, consider for example an ordinary gas in a spherical container of 
radius L, in which the number of composite particles is N, the mass of one particle is 
m and the mean velocity of particles is v. The ordinary statistical mechanics, with- 
out the help of thermodynamics, yields the partition function Zgas = Zgas{L, N, m, v) 
as simply a function of "parameters", L, N , m and v. Statistical mechanics, solely, 
can not determine what combinations of those parameters behave as state variables. 
To determine it, the first law of thermodynamics is necessary. (Note that the notion 
of heat in the first law is established by purely the argument in thermodynamics, 
not in statistical mechanics.) Comparing the differential of partition function with 
the first law results in the identification of partition function with the free energy 
divided by temperature. Then, since the free energy of ordinary gases is a function 
of the temperature and volume due to the "thermodynamic" argument, the partition 
function Zgas{L, N, m, v) should be rearranged to be a function of temperature and 
volume Zgssiy^T), where V = (47r/3)L^ and T = mv^ for ideal gases due to the law 
of equipartition of energy *\ (The dependence on N is, for example, Zgas oc N for 
ideal gases.) 

Here note that the reason why the temperature and volume are regarded as the 
state variables of the gas is that they are consistent with the four laws of "thermo- 



*' When the number of particles iV changes by, for example, a chemical reaction and an exchange 
of particles with environment, N is also the state variable on which the free energy depends. 



6 



Hiromi Saida 



dynamics" and have the appropriate properties as state variable. The appropriate 
properties are that the state variables are macroscopically measurable, the state vari- 
ables arc classified into two categories, intensive variables and extensive variables, 
and the extensive variables are additive. Those properties of state variables are 
specified by purely the argument in thermodynamics, not in statistical mechanics. 
Therefore, from the above, it is recognized that statistical mechanics can not yield 
the partition function as a "function of appropriate state variables" without the help 
of thermodynamics which specifies the appropriate state variables for the partition 
function. 

Turn our discussion to the Euclidean action method for curved spacetimes. Since 
the Euclidean action method is the technique to obtain the "partition function" of 
the spacetime under consideration (sec Appendix A), it is necessary to specify the 
state variables before calculating the Euclidean action. In this section, referring to 
Schwarzschild and de Sitter canonical ensembles,^^'^) we introduce the minimal set of 
assumptions for SdS thermodynamics, which specify the appropriate state variables 
for the partition function. Also, the special role of cosmological constant is clarified, 
which is already found in previous works.^^' The calculation of Euclidean action 
is carried out not in this section but in next section. 

2.2. SdS spacetime 

Before introducing the minimal set of assumptions, let us summarize the Lorentzian 
SdS spacetime in order to prepare some quantities used in the following discussions. 
The metric of SdS spacetime in the static chart is 

ds-' = -fir)dt'' + ^^+r''dn\ (2-1) 

fir) 

where df2^ = dO^ + sin^ 9 dip"^ is the line clement on the unit two-sphere, and 

f^r):=l-^-H'r^ , 3H^:=A, (2-2) 

where M is the mass parameter of black hole and A is the cosmological constant. 
The Penrose diagram of SdS spacetime is shown in Fig.l, and the static chart covers 
the region I. 

An algebraic equation /(r) = has one negative root and two positive roots. 
The smaller and larger positive roots are, respectively, the radius of BEH rj, and that 
of CEH Vc- The notion of CEH is observer dependent and the CEH at Vc is associated 
with the observer going towards the temporal future infinity i"*" in region I.-*^*^) The 
equation /(r) = is rearranged to 4f^ — 3f + s/O^ M H = 0, where f := ■\/3Hr/2. 
Then via a formula, sin^ = —4 sin^(^/3) -|- 3 sin(6'/3), we get 

2 /a\ 2 . /a + 2TT\ 
rh = —p^ — sm — , Tc = —p^ — sm , (2-3) 



where a is defined by, sin a := \/27 M H. The existence condition of BEH and CEH 
is < ^/27 M H < 1. This is equivalent to, < a < 7r/2, which means 

2M<n<3M< < < ;^ • (2-4) 
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Fig. 1. Penrose diagram of SdS spacetime. are the future/past null infinity, i'^ is the future 
temporal infinity. BEH is at r = rj, and CEH at r = r-c. Spacetime singularity is at r = 0. Static 
chart covers the region I. Semi-global black hole chart covers the regions I, 11;,, Illb and IV^. 
Semi-global cosmological chart covers the regions 1, flc, Illc and IVc. The maximal extension 
is obtained by connecting the two semi-global charts alternately. 



This denotes that r;, is larger than the Schwarzschild radius 2A4 and Tc is smaller 
than the de Sitter's CEH radius H~^. 

SdS spacetime has a timelike Killing vector := N dt, where is a normalization 
constant. This ^ becomes null at BEH and CEH. This means those horizons are 
the Killing horizons of S,. The surface gravity of BEH Kf, and that of CEH Kc are 
defined by the equations, V^^^\r=r,, = Kb^f^\r=r^ , V^^'^\r=rc = Kc^'^\r=rc- The 
surface gravity depends on N. Throughout this paper we take the normalization 
= 1 for BEH, and = — 1 for CEH to make Kc positive *-*. Then the surface 
gravities become equal to the absolute value |(l/2) df{r)/dr\ at each Killing horizon, 

^b = ^ (rc-n) (2rfe + re) = — {l-SH^l) , 

2n 2rb (2.5) 

Kc = — (rc-rfe) {rt + 2r,) = — (SH^r^^-l) , 

where Eq.(2-3) is used in the second equality in each equation. From the inequality 
ri) < rc in Eq.(2-4), we get 

Kb> Kc. (2-6) 

This implies the Hawking temperature of BEH is higher than that of CEH, which 
will be verified by Eqs.(4T) and (5T). 

For later use, let us show some differentials, 

l\/r fir^. 

(2-7a) 
(2.7b) 



drb _ 


2 


drb 


H 


M drb 


dM 


l-3H^rl 


dH ~ 


H dM 


Ore _ 


2 


Ore 


— - + 

H 


M drc 


dM 


3 r2 - 1 ' 


dH ~ 


H dM 



and 



dKb _ 1 l + Sff^rg dKb _ l-3gVg M dKb 

dM rll-^H'^rl ' dH 2Hrb H dM ' ^ ' ^' 



Even if TV = 1 for CEH, we can keep consistency of our analysis by changing the signature of 



Kc appropriately. 
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dKc _ 1 S-g^rg + l Okc 
'dM ~ ~r^3Wrf^ ' 'dW 

c 



2Hrr 



+ -T7 



H dM ' 



(2-7d) 



where we used a formula, cos 9 = 4 cos^(^/3) — 3 cos(^/3), and the differentials, 
dMC( = VTIH/ cos a, and dnot = ^fTJMj cos a, obtained from the definition of a, 
sin a := V^MH. 

The metric in semi-global black hole chart is given by the coordinate transfor- 
mation from {t,r,6,ip) to {rjbjXbjG^V^)'- 



where dr* := dr/f{r) which means 
2r* = In 



Vb + Xb :=e«''(*+'^* 





I /Kb 




I /Kb 






-In 




+ ln 












n + Tc 



l/Ka-l/Kb 



We get by this transformation, 

ds"" = nir) [ -dvl + dxl]+ dQ^ , 

where 



2M 

2 

Ktr 



r 

1 



l + Kb/Kc 



n + Tc 



1 



2-K(,/kc 



(2-8) 

(2-9) 

(2-10) 
(2-11) 



The transformation (2-8) implies the range of coordinates, —Xb < Vb < Xb a-^d < Xb) 
which covers the region I in Fig.l. By extending to the range, — oo < rfb < oo and 
— oo < X6 < oo, the semi-global black hole chart covers the regions I, lib, Hit and 
lYb in Fig.l. In these regions we find > since r < Tc- 

The metric in semi-global cosmological chart is given by the coordinate trans- 
formation from {t,r,9,(p) to {r]f.-,Xci^^v)'- 



,Kc(t-r*) 



% - Xc := e--^^" ■ ' , r]c + Xc-=-e 
where r* is given in Eq.(2-9). By this transformation we get 

ds'^ = Tcir) [ -dril + dxl]+ dQ"^ , 

where 



Kc{t+r*) 



2M 



K^r \rb 



l + Kc/Kb 



n + Tc 



+ 1 



2-Kc/Ki, 



(2-12) 

(2-13) 
(2-14) 



The transformation (2-12) implies the range of coordinates, Xc < Vc < ~Xc and 
Xc < 0, which covers the region I in Fig.l. By extending to the range, — oo < r)c < oo 
and — oo < Xc < oo, the semi-global cosmological chart covers the regions I, lie, IHc 
and IV c in Fig.l. In these regions we find Tc> since rb < r. 

The maximally extended SdS spacetime is obtained by connecting the two semi- 
global charts alternately as shown in Fig.l. 
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2.3. Minimal set of assumptions and working hypothesis 

As mentioned in Sec. 2.1, we introduce the minimal set of assumptions with re- 
ferring to Schwarzschild thermodynamics formulated by York.^^ Those assumptions 
should construct thermal systems for SdS thermodynamics and give us enough state 
variables which appear as independent variables in the free energy of our thermal 
system. There are three key points in Schwarzschild thermodynamics from which we 
can learn about the way to ensure the "thermodynamic consistency" in SdS ther- 
modynamics. Those key points are the same from which the basic assumptions of 
de Sitter thermodynamics are introduced in previous paper, and we summarize 
those three key points in Appendix B of this paper. 

Here we must comment on Anti-de Sitter (AdS) black holes.^^ AdS black hole 
thermodynamics has a conceptual difference from the other black hole thermodynam- 
ics. The difference appears, for example, in the definition of temperature. While the 
temperatures in asymptotic flat black hole and de Sitter thermodynamics include 
the Tolman factor the temperature assigned to AdS black hole^^ does not 
include the Tolman factor, where the Tolman factor^^) expresses the gravitational 
redshift affecting the Hawking radiation propagating from horizon to observer (see 
for example Eq.(4-1) in Sec.4.1 and the key point 3 of Schwarzschild thermodynamics 
in Appendix B). The temperature in AdS black hole thermodynamics can not be 
measured by a thermometer of the physical observer who arc outside the black hole. 
This implies that the state variables in AdS black hole thermodynamics are defined 
not by the observer outside the black hole, but defined just on the black hole event 
horizon on which no physical observer can rest. In this paper we do not refer to AdS 
black hole thermodynamics, since it seems to be preferable to expect that state vari- 
ables are defined according to a physical observer. Hence we refer to Schwarzschild 
thermodynamics, which is based on a physical observer (See Appendix B). 

2.3.1. Zeroth law and independent variables 

We will construct two thermal equilibrium systems for BEH and CEH, but place 
only one observer who can measure the state variables of BEH and CEH. Such 
observer is in the region 1, rh < r < rc (see Fig.l). However, as mentioned at 
Eq. (2-6) , Hawking temperature of BEH is higher than that of CEH. This temperature 
difference implies that, when the region I constitutes one thermodynamic system, 
the thermodynamic state of region I is in a non-equilibrium state. Therefore, by 
dividing the region I into two regions, we construct two thermal equilibrium systems 
for BEH and CEH individually which are measured by the same observer. To do so, 
we adopt the following assumption as the zeroth law: 

Assumption 1 (Zeroth law) Two thermal "equilibrium" systems for BEH and 
CEH in SdS spacetime are constructed by the following three steps: 
1. Place a spherically symmetric thin wall at r = in the region I, r^, < < 
as shown in Fig. 2. This wall has negligible mass, and reflects perfectly Hawking 
radiation coming from each horizon. We call this wall the "heat wall" hereafter. 
The BEH side of heat wall is regarded as a "heat bath" of Hawking temperature 
of BEH due to the reflected Hawking radiation. Also the CEH side of heat wall 
is regarded as a heat bath of Hawking temperature of CEH. 
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Dh:rh<r< r,, /- - : r,, <r<rc 
1+ 




Heat Wall : r = Vw (Observer is here) 



Fig. 2. Our thermal equilibrium systems for BEH (-Dt) and CEH (-De)- State variables of them are 
defined at the heat wall. 



2. The region enclosed by BEH and heat wall, ri, < r < r^, is filled with 
Hawking radiation emitted by BEH and reflected by heat wall, and forms a 
thermal equilibrium system for BEH. Similarly the region Dc enclosed by CEH 
and heat wall, < r < rc, forms a thermal equilibrium system for CEH. 
Hence we have "two" thermal equilibrium systems separated by the heat wall. 
In the statistical mechanical sense, these two thermal equilibrium systems are 
described by the canonical ensemble, since those systems have a contact with 
the heat wall. 

3. Set our observer at the heat wall. When the observer is at the BEH side of heat 
wall, the observer can measure the state variables of thermal equilibrium system 
for BEH. The same is true of CEH. Then the state variables of two thermal 
equilibrium systems are defined by the quantities measured by the observer at 
heat wall. 

Note that the two thermal equilibrium systems constructed in this assumption have 
ah'eady been used by Gibbons and Hawking^*^) to calculate Hawking temperatures 
of BEH and CEH. Also the above step 3, which gives a criterion of defining state 
variables, has already been adopted in the consistent thermodynamics of single- 
horizon spacetimes.^)"^-* We can regard this assumption as a simple extension of the 
key point 1 of Schwarzschild thermodynamics shown in Appendix B. 

It is expected that state variables of the thermal equilibrium system for BEH 
depend on BEH radius rf, and/or BEH surface gravity Kh- Similarly, state variables 
of CEH depend on rc and/or Kc. These imply that the state variables of BEH and 
CEH depend on M and A, since the horizon radii and surface gravities depend on 
M and A via Eqs.(2-3) and (2-5). Furthermore, by the step 3 in assumption 1, there 
should be r^-dependence in state variables of BEH and CEH, since the observer is at 
r = r^j. Therefore the state variables of BEH and CEH depend on three parameters 
M, A and r^. 

The existence of three parameters (M, yl, r^) may imply that the CEH is re- 
garded as a source of external gravitational field which affects the thermodynamic 
state of BEH. Also, BEH is a source of external gravitational field affecting the ther- 
modynamic state of CEH. Here it is instructive to compare qualitatively our two 
thermal equilibrium systems of horizon with a magnetized gas. The gas consists of 
molecules possessing a magnetic moment, and its thermodynamic state is affected 
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by an external magnetic field. The qualitative correspondence between the magne- 
tized gas and our thermal equilibrium systems of horizon is described as follows; the 
gas corresponds to the system Df, (Dc), and the external magnetic field corresponds 
to the external gravitational field produced by CEH (BEH). Then, what we should 
emphasize is the following fact of the magnetized gas: When the gas is enclosed in 
a container of volume V^as and an external magnetic field i?ex is acting on the gas, 
the free energy Fg^s of the gas is expressed as a function of three independent state 
variables (see for example §52, 59 and 60 in Landau and Lifshitz^^^), 

-^gas ~ -^gas(^gas) ^as) -f^ex) ; 

where Tgas is the temperature of the gas. According to this fact of the magnetized 
gas, it is reasonable for our two thermal equilibrium systems of horizons to require 
that the free energies are also functions of three independent state variables. For the 
BEH, the free energy Fi, is 

Fb = Fb{Tb,Ab,Xi,), (2-16) 

where Tf, is the temperature of BEH, is the state variable of system size, and 
is the state variable which represents the effect of CEH's gravity on the BEH. And 
for the CEH, the free energy Fc is 

Fc = F,{n,A„X,), (2-17) 

where Tc is the temperature of CEH, Ac is the state variable of system size and 
Xc is the state variable which represents the effect of BEH's gravity on the CEH. 
Indeed, it will be proven in Sees. 4.1 and 5.1 that the thermodynamic consistency 
never hold unless the free energies are functions of three independent variables as 
shown in Eqs.(2-16) and (2-17). 

Now we recognize that, because free energies are functions of three independent 
variables (as will be verified in Sees. 4.1 and 5.1), the following working hypothesis 
is needed: 

Working Hypothesis 1 (Three independent variables) To ensure the thermo- 
dynamic consistency of our thermal equilibrium systems constructed in assumption 1, 
we have to regard the cosmological constant A as an independent working variable. 
Then the three quantities {M^A^r^) are regarded as independent variables, and con- 
sequently the free energies Ff, and Fc of our thermal equilibrium systems are functions 
of three independent state variables as shown in Eqs. (2-16) and (2T7). On the other 
hand, when we regard the non-variable A as the physical situation, it is obtained by 
the "constant A process " in the consistent thermodynamics for BEH and CEH which 
are constructed with regarding A as a working variable. 

This working hypothesis will be verified in Sees. 4.1 and 5.1, and we can not 
preserve thermodynamic consistency without regarding A as an independent working 
variable. This implies that, as already commented in Sec.V of previous paper, 
the thermal equilibrium states of event horizon with positive A may construct the 
"generalized" thermodynamics in which A behaves as a working variable and the 
physical process is described by the constant A process. 



(2-15) 



12 



Hiromi Saida 



2.3.2. Scaling law and system size 

In thermodynamics of ordinary laboratory systems, all state variables are clas- 
sified into two categories, extensive state variables and intensive ones. The criterion 
of this classification is the scaling behavior of state variables under the scaling of sys- 
tem size. However, as explained by the key point 2 of Schwarzschild thermodynamics 
shown in Appendix B, the state variables in thermodynamics of single-horizon space- 
times^^"^) have its own peculiar scaling behavior classified into three categories, and 
the state variable of system size is not a volume but the surface area of heat bath. 
Although the scaling behavior differs from that in thermodynamics of ordinary lab- 
oratory systems, the peculiar scaling behavior in thermodynamics of single-horizon 
spacetimes retains the thermodynamic consistency as explained in the key point 3 
in Appendix B. Then, we assume that the key point 2 of Schwarzschild thermody- 
namics is simply extended to our two thermal equilibrium systems constructed in 
the assumption 1: 

Assumption 2 (Scaling law and system size) AH state variables of our ther- 
mal equilibrium systems are classified into three categories; extensive variables, in- 
tensive variables and thermodynamic functions, and satisfy the following scaling law: 
When the length size L (e.g. horizon radius) is scaled as L ^ XL with an arbitrary 
scaling rate A (> 0), then the extensive variables X (e.g. system size) and intensive 
variables Y (e.g. tem,perature) are scaled respectively as X ^ /(^ X and Y Y , 

while the thermodynamic functions 'P (e.g. free energy) are scaled as ^ 
This implies that the thermodynamic system size of our thermal equilibrium systems 
should have the areal dimension, since the system size is extensive in thermodynamic 
argument. Then we assume that the surface area of heat wall, A := 47rr^, behaves as 
the consistent extensive variable of system size for our thermal equilibrium systems 
for BEH and CEH. This denotes to set Ab = Ac = Ain Eqs.{2-16) and (2-17). 

Accepting this assumption, the length size scaling in our thermal equilibrium 
systems for BEH and CEH should be specified. Here recall that, due to the working 
hypothesis 1, the fundamental independent variables in our thermal equilibrium 
systems are M, A and ryj. Therefore the fundamental length size scaling for our 
thermal equilibrium systems is composed of the following three scalings; 

M^XM , H^^H , r^^Xryj, (2-18) 
A 

where H is defined by 3H^ := A in Eq.(2-2), and A is an arbitrary scaling rate. The 
extensivity and intensivity of each state variable of our thermal equilibrium systems 
should be defined under these fundamental length size scalings as explained in the 
assumption 2. *) 



*^ As explained in Appendix B in previous paper,'"' when we regard ^ := 4 7rr^ as a state 
variable of system size, the scaling of system size should be restricted to the homothetic one, which 
is the spherical scaling due to the spherical symmetry of SdS spacetime. The fundamental length 
size scaling (2-18) is consistent with this restriction. See Appendix B in previous paper''' for details 
of such restriction. 
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2.3.3. Euclidean action method and how to obtain state variables 

We need to specify how to get the state variables. As noted in the step 2 
in assumption 1, thermodynamics of our two thermal equilibrium systems should 
be constructed in the canonical ensemble. Therefore we use the Euclidean action 
method which is the technique to calculate the partition function of quantum grav- 
ity.^^^ Indeed, the Euclidean action method has already made successes to obtain the 
partition function of canonical ensemble for the thermodynamics of single-horizon 
spacetimes.'^^"'^^ The key point 3 in Appendix B explains how to use the Euclidean 
action and to define state variables for Schwarzschild thermodynamics. Then, refer- 
ring to the key point 3, we adopt the following assumption: 

Assumption 3 (State variables and Euclidean action method) Euclidean ac- 
tions I^b o-i^d Iec of our two thermal equilibrium, systems yield the partition functions 
of canonical ensembles by Eq. (A-3) of Appendix A. And the free energies F}, and Fc 
are dejined by Eq. {A-5), where the temperatures are defined by i?g.(A-4). Then, once 
Ff, and Fc are determined, all state variables of BEH and CEH are defined from Ff, 
and Fc as for thermodynamics of ordinary laboratory systems. For example, BEH 
entropy Sb is defined by Sb ■= —dFb/dTb, where Tb is the temperature of BEH. 

As seen in Eq.(A-2), the Euclidean action is produced from the Lorentzian action. 
The Lorentzian Einstein-Hilbert action II is 



where M is the spacctimc region under consideration, R is the Ricci scalar, g is the 
determinant of metric, h and K in the second term are respectively the determinant 
of first fundamental form (induced metric) and the trace of second fundamental 
form (extrinsic curvature) of the boundary dAi, and Iq is the integration constant 
of II and called the subtraction term. The second term in Eq.(2-19) is required to 
eliminate the second derivatives of metric from the action. ^^-^ The Iq is independent 
of the metric in M, and does not contribute to the Einstein equation obtained by 
SIl = 0. The Einstein equation for SdS spacctime gives the relation R = 4A. 

When we use the Euclidean action method, it is necessary to specify the inte- 
gration constant Iq. It is natural to require that our thermal equilibrium systems for 
BEH and CEH should reproduce, respectively, the Schwarzschild thermodynamics 
in the limit A ^ and the de Sitter thermodynamics in the limit M — 0. Then the 
following working hypothesis is naturally required: 

Working Hypothesis 2 (Integration constants in Euclidean action) For the 

thermal equilibrium system for BEH, the integration constant in Euclidean action is 

determined with referring to Schwarzschild canonical ensemble formulated by York.^^ 
For the thermal equilibrium system, for CEH. the integration constant in Euclidean 
action is determined with referring to de Sitter canonical ensemble formulated in 
previous paper. 

We introduce this working hypothesis as if this is a statement separated from 
the assumption 3. However the determination of integration constant accompanies 
necessarily the Euclidean action method. The working hypothesis 2 is regarded as a 
part of the assumption 3. 




(2-19) 
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2.3.4. Effects of external gravitational fields 

By the assumption 3 together with the working hypothesis 2, the concrete func- 
tional form of free energies Fi, and Fc can be determined as functions of three inde- 
pendent working variables (M, yl, r^). However, since the form of the state variables 
Xh and Xc have not been specified yet, we can not rearrange Ff, and Fc to functions 
of independent state variables, {Th,A,Xf)) and {Tc, A, Xc). 

As explained at Eqs.(2-16) and (2-17), the CEH (BEH) is regarded as the source 
of external gravitational field which affects the thermodynamic state of BEH (CEH). 
This means that the state variables X^ and Xc represent, respectively, the thermo- 
dynamic effect of CEH's gravity on BEH and that of BEH's gravity on CEH. Then 
it is reasonable to expect that X{, depends on the quantity characterizing the CEH's 
gravity, and Xc depends on the quantity characterizing the BEH's gravity. Moreover, 
due to the step 3 in assumption 1, Xf, and Xc should be measurable for the observer 
at rw Then, we can offer two candidates for the pair of dimensionless characteristic 
quantities of BEH's and CEH's gravities; 

• First candidate pair consists of Kf,r^ and Kcr^y, where Kf,r^ is for BEH's gravity 
and KcVyj is for CEH's gravity. This pair means that both of BEH's and CEH's 
gravities are characterized by three quantities (M, yl, r^), since Kf, and Kc depend 
on M and A. 

• Second candidate pair consists of M/r^j and Hr^^, where M/r^j is for BEH's 
gravity and Hvy, is for CEH's gravity. This pair means that the BEH's gravity 
is characterized by (M, r^), and the CEH's gravity is characterized by {A,rw). 

Here, purely logically, we can consider the "inverse" pair of second one, where 
is for BEH's gravity and M/r^ is for CEH's gravity. This means that the BEH's 
gravity is characterized by {A, Tw), and the CEH's gravity is characterized by (M, r^). 
However this is physically unacceptable, since we do not expect that the BEH does 
not depend on M and the CEH docs not depend on A. Therefore, the reasonable 
candidates for the pair of characteristic quantities of BEH's and CEH's gravities 
are the two candidates listed above. Then, should be a function of (kc , r^) or 
{H , ryj), and Xc should be a function of , r^) or (M , ryj). 

On the other hand, as will be mathematically verified in Sees. 4. 3 and 5.3, the 
state variables Xi, and Xc are the extensive variables and proportional to r^. The 
proportionality to is consistent with the scaling law of extensive variables denoted 
in assumption 2. And, according to the previous paragraph, the factor of proportion- 
ality should be a function of the characteristic quantity of BEH's or CEH's gravity. 
Although the verification of the extensivity of Xf, and Xc are shown later, we accept 
it in the following assumption 4 for the simplicity of our discussion. 

From the above, it is reasonable to adopt the following assumption: 
Assumption 4 (Extensive variable of "external field" ) The state variables X^ 
and Xc in Eqs. {2-16) and (2-17) are the extensive variables. (This will be verified 
in Secs.A.S and 5.3). Then, there are two candidates for the functional forms of X^ 
and Xc. One of them is based on the quantities {n},rw , Kcry,): 

Xb = rl^bi^cryj) , Xc = rl^ci^bryj) , (2-20) 

where •P'b and <Pc are arbitrary functions of single argument, whose functional forms 
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are not specified at present. Another candidate of Xf, and Xc is based on the quan- 
tities {M/r-uj , Hr-u,): 

Xi, = rl%{Hr^) , X, = rl^,{M/r^). (2-21) 

At least for the present author, no criterion to choose one of these candidates is 
found, and the way for determining the functional forms of and cire also un- 
known. An obvious constraint on and Ifc is that they never be constant to make 
Xf) and Xc independent of the sta,te variable of system size A := 47rr^. 

In Sec. 6.1, we will make some comments on the issue which of Eqs.(2-20) and (2-21) 
is valid. Those comments will suggest that Eq.(2-20) may be appropriate, but we 
do not have mathematical verification to choose Eq.(2-20) as the general form of 
Xb and Xc- Therefore, to retain the logical strictness of this paper, we list the two 
possibilities (2-20) and (2-21) in the assumption 4. 

Prom the above, we recognize that the minimal set of assumptions for "consis- 
tent" thermodynamics of our two thermal equilibrium systems should be composed 
of four assumptions. However, the determination of functions and remains 
as a future task and we can not find concrete functional forms of them. Although 
the state variables Xj, and Xc are not specified in this paper, the existence of them 
enables us to examine the validity of entropy-area law in SdS spacetime as shown in 
Sees. 4 and 5. 



§3. Euclidean actions 



Referring to the assumption 3 and working hypothesis 2, we calculate Euclidean 
actions for the two thermal equilibrium systems for BEH and CEH constructed in 
the assumption 1. 

3.1. Euclidean action for BEH 

Euclidean space of thermal equilibrium system for BEH is obtained by the Wick 

rotations t —ir in the static chart and r/5 — > —iLUb in the semi-global black 
hole chart. These Wick rotations arc equivalent, because the coordinate transforma- 
tion (2-8), r]b = e*^^^ sinh (k^ i), implies that the imaginary time in the semi-global 
chart is defined by w;, := e^^'^* sin (/tftr), where r is the imaginary time in the static 
chart. Euclidean metric in the static chart is 

ffv^ 

dsl = fir)dT^ + j^ + r^df2\ (3-1) 

Euclidean metric in the semi-global black hole chart is 

dsl = n{r) [ dLjl + dxb] + dn'' , (3-2) 

where is defined in Eq.(2-ll). About the semi-global chart, we get from the 
coordinate transformation (2-8), 
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Then, because the thermal equilibrium system for BEH is the region Df,, < r < r^, 
in Lorentzian SdS spacetime, we find the topology of Euclidean space of thermal 
equilibrium system for BEH is x S'^. Because o;^ + = at r = r^, the center 
of i:'2.part is at the BEH r = ri,, and the boundary of D^-part is at the heat wall 
r = Tyj. The topology of heat wall boundary is x S"^, where is along the 
T-direction. 

Because the Lorentzian SdS spacetime is regular at r^, the Euclidean space is 
also regular at r;,. To examine the regularity of Euclidean space at r;,, we make use 
of the static chart (3T). Let us define a coordinate and a function ^{x^) by 

xl:=r-n , l{xb) ■■= \Jf{n + xl) . (3-4) 
We get ^'{xb) ■= d'y{xb)/dxb = [xb/^{xb)] df{r)/dr, from which we find 



1- u ^ 



xb^o 7(xb) hm 7 [xb) 



This means 7'(0) = ^/2l^, and near the BEH, / ~ [7(0) + 7'(0) xf = 2 Kxg. There- 
fore the Euclidean metric near BEH is 

ds%^— \ xl d{Kb rf + dxl]+ df2^ . (3-6) 

Kb 

It is obvious that the Euclidean space is regular at BEH if the imaginary time has 
the period defined by 

< T < /?6 := — . (3-7) 

Kb 

Throughout our discussion, r has the period Pb in the Euclidean space of thermal 
equilibrium system for BEH. 

Now let us proceed to the calculation of the Euclidean action lEb of thermal 
equilibrium system for BEH. Following the working hypothesis 2, wc use the same 
integration constant Iq as Schwarzschild canonical ensemble,^)' which gives us 

/o := -/f-t) = -± / dx^VhK^^-'\ (3-8) 
oTT Jqm 

where 7^^^*^ is the Lorentzian Einstein-Hilbert action for fiat spacetime which reduces 

to only the surface term (the second term in Eq.(2-19) ) due to i? = and yl = 
for flat spacetime, and i^^''^*) is the trace of second fundamental form of dAA in flat 
spacetime. ^^'^^^ Here note that, the integral element y/h in J^^^*^ should be given 

by that of SdS spacetime when is used as the integration constant of action 

integral of SdS spacetime, because the background spacetime on which the integral 
in J^'''^*) is calculated is SdS spacetime. 

For our thermal equilibrium system for BEH in SdS spacetime, the region M in 
II is Db, Tb < r < Tw, and its boundary dDb is at r^j. There is another boundary at 
rb in the Lorentzian region Db- However we do not need to consider it, because the 
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points at rj, in the Euclidean space do not form a boundary but are the regular points 
when r has the period (3-7). Then the first fundamental form hij {i,j = 0,2,3) of 
dDi, in the static chart is 



hij dx^ dx^ = -fw dt^ + rl, dQ"^ 



where 



fw ■■= f{rw) = 1 



2M 



(3-9) 



(3-10) 



Here, since D^, is the region enclosed by BEH and heat wall, the direction of unit 
normal vector to dDi, is pointing towards CEH, n'* oc dr- Then the second 
fundamental form of dDj, in the static chart is 



\A«diag. 



— TT + H^r. 



W J ' W 1 ' w 



sin 



(3-11) 



where diag. means the diagonal matrix form. 

On the other hand, the second fundamental form K^^^^ of a spherically sym- 
metric timelike hypcrsurface of radius in flat spacetime is given by setting M = 
and = in Eq.(3-ll), 



^(flat) ^ [Q,ryj,ryj so? e ] 



(3-12) 



This gives i^^fl'^*) := K'^ K^f"^'^ = 2r-^ 

Prom the above, applying the Wick rotation t — ^ —ir to the Lorentzian action II 
in Eq.(2-19), we obtain the Euclidean action of the thermal equilibrium system 
for BEH via Eq.(A-2), 



3H^ 
"StT 

2 



dx'^ 



'9E + ^ [ 

OTT Jq 



dx% 



DEb 



IdDEb 



hE {^Ke 



K 



(flat) 



(3-13) 



where the relation for SdS spacetime R = AA = 12 is used in the first equality, 
the relation M — H^rf = 3M — Vb due to /(r^) = is used in the second equality, 
Qe the quantity Q evaluated on Euclidean space, and DEb is the Euclidean region 
denoted by < r < /?5 , < r < r^t,. , < < vr and < ip < 2ir. This 
lEb corresponds to Ie[9Eci\ in Eq.(A-5) of Appendix A, which yields the partition 
function of our thermal equilibrium system for BEH. 

Note that lEb should reproduce the Euclidean action of Schwarzschild canonical 
ensemble as required in the working hypothesis 2. To check if this is satisfied, take 
the limit A ^ 0, 



lim lEb = 47rM 



M- 2r„ 



A=0 



(3-14) 



This coincides with the Euclidean action of Schwarzschild canonical ensemble for- 
mulated by York.^^ 
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3.2. Euclidean action for CEH 

Euclidean space of thermal equilibrium system for CEH is obtained by the Wick 
rotations t ^ —it m. the static chart and rjc — ^ —i <^c in the semi-global cosmological 
chart. These Wick rotations are equivalent, because the coordinate transforma- 
tion (2-12), r/c = 6^"'=''* sinh(Kct), implies that the imaginary time ojc in the semi- 
global chart is defined by u)c '■= e""''^ sin(Kcr), where r is the imaginary time in 
the static chart. Euclidean metric in the static chart is given by Eq.(3-1). Euclidean 
metric in the semi-global cosmological chart is ds^ = Tc(r) [dcj^ + ^Xc] +f'^dQ'^^ 
where Tc is defined in Eq.(2-14). About the semi-global chart, we get from the 
coordinate transformation (2-12), 

Then, because the thermal equilibrium system for CEH is the region Dc, rw < r < Tc, 
in Lorentzian SdS spacetime, we find the topology of Euclidean space of thermal 
equilibrium system for CEH is x S^. Because o;^ + Xc = *t r = Tc, the center 
of -D^-part is at the CEH r = r^, and the boundary of D^-part is at the heat wall 
r = Tyj. The topology of heat wall boundary is x S'^, where is along the 
r-direction. 

Because the Lorentzian SdS spacetime is regular at Tc, the Euclidean space is 
also regular at Vc- Using the static chart (3T) and defining a coordinate and a 
function 7(xc) by x^. := rc — r and 'y{xc) '■= ^/ /{vc — x'^), we obtain the Euclidean 
metric near CEH, 

ds\c:±— [xld{KcTf + dxl]+r'^dfi'^ . (3-16) 

It is obvious that the Euclidean space is regular at CEH if the imaginary time has 
the period /3c defined by 

0<T</3c:=— . (3-17) 

Throughout our discussion, r has the period j3c in the Euclidean space of thermal 
equilibrium system for CEH. 

Now we proceed to the calculation of Euclidean action Iec of the thermal equi- 
librium system for CEH. Lorentzian action is defined in Eq.(2-19) for the spacetime 
region Dc, < r < rc- The calculation of Euclidean action I^c is parallel to that 
of lEb except for the direction of unit normal vector to dDc and the integration 
constant Iq in 1^. Concerning the vector n'^, since Dc is the region enclosed by CEH 
and heat wall, the direction of is pointing towards BEH, oc —dr- 

Concerning the integration constant, following the working hypothesis 2, we de- 
termine the integration constant /q for CEH with referring to the de Sitter canonical 
ensemble formulated in previous paper.^) At the limit M —>■ 0, the term Jq should 
reduce to the integration constant /q^^^ of the de Sitter canonical ensemble. 



Io(M 



0) = 1^ := (^j^ - l) V^(M = 0)/f"*) , (3-18) 
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where Ii^^ is the action of flat spacetime used in Eq.(3-8). Note that the CEH 
radius in de Sitter spacetime is H~^, and the factor H r^ is the ratio of heat wall 
radius ryj to CEH radius. Hence we set Iq for the CEH in SdS spacetime, 

lo := - l) v^/f^) , (3-19) 

where it should be emphasized that the signature of K^^^^^ in the integrand of J^^^^^ 
(shown in Eq.(3-8)) should be reversed, because the direction of normal vector is 
reversed as mentioned in previous paragraph. 

Then we obtain the Euclidean action Iec of our thermal equilibrium system for 
CEH via Eqs.(2-19) and (A-2), 

^£;c = -Y [3M-re + 2rc/^] , (3-20) 

where the relation M—H^r^ = SM—Vc due to f{rc) = is used, and the overall minus 
signature comes from the direction of normal vector n'^ to dDc- This Iec corresponds 
to lEidEd] in Eq.(A-5), which yields the partition function of our thermal equilibrium 
system for CEH. 

Note that Iec should reproduce the Euclidean action of de Sitter canonical en- 
semble as required in the working hypothesis 2. To check if this is satisfied, take the 
limit M ^ 0, 

hm /Ee = -;^ [l-2iHr^f] . (3-21) 

This coincides with the Euclidean action of de Sitter canonical ensemble formulated 
in previous paper. ^) 



§4. Black hole event horizon 



We examine whether the entropy-area law holds for "consistent" thermodynam- 
ics of our thermal equilibrium system for BEH. 

4.1. Temperature and free energy of BEH 

By the assumption 3, the temperature Tf, of BEH is defined by Eq.(A-4) of 
Appendix A, which relates Ti, to the proper length in the imaginary time direction 
at the boundary dDj, (the direction along part of boundary topology x in 
Euclidean space). 



Tb:= 



-, -1 



, dr 



Kb 



27r 



(4-1) 



where Pb is the imaginary time period (3-7) and fw is in Eq.(3-10). Under the length 
size scaling (2-18), this temperature is scaled as Tb ^ X~^Tb. Therefore, by the 
assumption 2, T5 is an intensive state variable of BEH. 

Note that this coincides with the Hawking temperature of BEH derived orig- 
inally by Gibbons and Hawking, ^'^^ and the factor ^/f^ is the so-called Tolman 
factor^'^) which expresses the gravitational redshift affecting the Hawking radiation 
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propagating from BEH to observer at r^. Therefore this Tj, is the temperature 
measured by the observer at heat wall. 

By the assumption 2, the extensive state variable of system size for our thermal 
equilibrium system is the surface area of heat wall, 



A := Airri 



(4-2) 



By the assumption 3, the free energy of BEH in Eq.(2-16) is defined by 
Eq.(A-5), 

F^{n, A, Xfe) := -n lEb = - VX - ^ ' ■ (4-3) 



2 v/u) 

Under the length size scaling (2-18), this free energy satisfies the scaling law of 
thermodynamic functions, Fh XFh. As discussed at Eq.(2-16), Ff, is regarded as 
a function of three independent state variables Tf,, A and the state variable of 
CEH's gravitational effect on BEH. However, since is not specified as mentioned 
in the assumption 4, the form of Fi, as a function of {Ti,, A, Xi,) remains unknown. 
Instead, Eq.(4-3) shows F^ as a function of independent parameters (ikf , H , r^y). 

Let us verify that the free energy Fj, is a function of three independent state 
variables. We use the reductive absurdity: Assume that only two (not three) 
state variables are independent. This assumption means that, as for the ordinary 
non-magnetized gases, F^ is a function of Tf, and A, Fiy{Tb,A). Here it is obvi- 
ous from Eq.(4-1) that depends on three parameters {M, H,ru)), while A de- 
pends only on r^y. Then, via Eq.(C-5) of Appendix C, a mathematical relation 
{dMFb)/{dMTb) = {dHFi,)/{dHTb) must hold if F^ is a function of {Tb,A). However 
we find this relation does not hold, {dMFh) / {dMTb) 7^ {duFb) / [dnTb)^ via Eq.(4-5) 
shown below. Hence the assumption of two independent state variables is denied by 
the reductive absurdity. Now the working hypothesis 1, which assumes F^ to be a 
function of three independent state variables, is verified. 

To support the discussion in previous paragraph and for later use, we show some 
differentials: 



dM 

dTb 



1 



27rr2/^/' 



n 



9H AirHrbf, 



3/2 

w 



(l-3i7V,^)-i + ^^'^^'^^^^ 



n 



and 



dFb 


1 




dM 


2/5' 




dFb 


1 


(2r^ +rb- 


dH 


" 2fJ' 



+ 



1 + 3H\^ 

Jii 



SH'^r 



where the differentials in Eqs.(2-7a) and (2-7c) are used. Then we get 

dFb 



dM 



2dTb 

= -TTTb 



dM 



(4-4a) 
(4-4b) 

(4-4c) 
(4-4d) 

(4- 5a) 



To what extent is the entropy- area law universal ? 



21 



where the definition of r^, /(r^) = 0, is used in the first relation. These relations are 
used in examining the entropy-area law in next subsection. 

4.2. Entropy of BEH 

By the assumption 3, the entropy S^, of BEH is defined as the thermodynamic 
conjugate variable to T^, 

dF,{T,,A,X,) 
S, .= . (4-6) 

Under the length size scaling (2-18), this entropy satisfies the extensive scaling law, 
Sb Sf,. By Eq.(C-3) of Appendix C, Sf, is rearranged to 

' (dMTb) (dnXb) - {dnn) (SmX^) ' ^' > 

Then we get by Eq.(4-5a), 

{dMn){dHXb)-{dHn){dMXb) ■ ^ 

This denotes the following: If djuXh = 0, then the entropy-area law 5*6 = vr holds. 
However, OmXi, ^ 0, then Eq.(4-5b) together with Eq.(4-8) imply that the entropy- 
area law breaks down Sb ^ T^r^. Therefore we find that the entropy-area law holds 
if and only if BmX^ = 0. 

In summary, although the validity of entropy-area law for BEH can not be judged 
at present, we can clarify the issue on the entropy-area law for BEH: The necessary 
and sufficient condition to ensure the entropy- area law for BEH is that the BEH is 
in thermal equilibrium and the state variable Xf, satisfies OmXi, = 0. // the CEH's 
gravitational effect on BEH, X^, is characterized by the CEH's quantity Kc^w as 
shown in Eq. {2-20), then OmX^ ^ and the entropy-area law breaks down for BEH. 
If the CEH's gravitational effect X^ is characterized by Hryj as shown in Eq.{2-21), 
then OmXi, = and the entropy-area law holds. The validity of entropy-area law for 
BEH will be judged by revealing which of Kc or H is appropriate as the characteristic 
quantity of CEH's gravity. 

4.3. Thermodynamic consistency of BEH 

The remaining part of Sec. 4 is for the "thermodynamic consistency" of our 
thermodynamics of BEH under the minimal set of assumptions. Since the concrete 
form of Xb is not specified, the following discussions and calculations are very formal. 
However we can imply that the thermodynamic consistency is satisfied. Moreover, 
it is also verified that Xb is an extensive variable and proportional to r^. 

Following the assumption 3, the internal energy Eb of BEH can be defined by 
the argument of statistical mechanics. 



Eb:=- 



A,Xi,=const. ^(V^b) 



^^P^l^^) =Fb + TbSb, (4-9) 
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where Eq.(A-5) of Appendix A is used in the second equahty, and the definition of Sb 
in Eq.(4-6) is used in the third equahty. Under the length size scahng (2-18), this Et, 
satisfies the scahng law of thermodynamic functions, Efy AE'b. The third equality 
in Eq.(4-9), Efy = + TbS{j, is regarded as the Legendre transformation between F^ 
and Eh, which determines Ef, to be a function of {Sb, A, Xh). On the other hand, in 
thermodynamic argument, the internal energy is the thermodynamic function which 
is regarded as a function of only extensive state variables. Hence, it is verified that 
Xb must be an extensive state variable. (The proportionality of Xb to will also 
be shown mathematically at the end of this subsection.) 

Next, in order to see the first law, we need the intensive state variables which 
are thermodynamically conjugate to A and Xb- By the assumption 3, the conjugate 
state variable is defined by the appropriate differential of a thermodynamic function. 
In an analogy with the ordinary pressure of ordinary gases, the surface pressure (at 
the heat wall) ab is defined formally as 



(Tb-= -- 



dFb{Tb,A,Xb) 



dA 



(4- 10a) 



1 



1 



(4- 10b) 



Svrr^ {duTb) (dnXb) - (^mX,) {dnn) 
{ {duFb) {dr^Tb) - (dr^Fb) {dnTb) } (duXb) 
+ { {dr^Fb) {duTb) - {duFb) {dr^Tb) } {dnXb) 
+ { {duFb) {dnTb) - (OffFb) {dun) } {drM 

where Eq.(C-2) of Appendix C is used in the second equality. Under the length 
size scaling (2-18), this ab satisfies the intensive scaling law, ab — \~^ab- The 
state variable given by the same definition with ab appears also in single-horizon 
thermodynamics^-*"^) to ensure the thermodynamic consistency. (See Appendix B in 
previous paper^) for thermodynamic meanings of A and ab-) 

The intensive state variable Yb conjugate to Xb is defined formally as 



Yb-r- 



dFb{Tb,A,Xb) _ jduFb) {dnTb) - jdsFb) {dun) 
dXb {dMXb) {dnTb) - {dnXb) (SmT;.) ' 



(4-11) 



where Eq.(C-3) of Appendix C is used in the second equality. Under the length size 
scaling (2-18), when Xb is scaled as an extensive variable Xb — ^ }? Xb, then this Yb 
satisfies the intensive scaling law, Yb A"^ Yj,. 
Then by definitions of Sb, ab and Yb, we get 

dFb{Tb, A, Xb) = -Sb dTb -abdA + Yb dXb . (4-12) 

The first law follows this relation via the Legendre transformation in Eq.(4-9), 

dEb{Sb, A, Xb) = Tb dSb -abdA + Yb dXb . (4-13) 



Concerning the internal energy, the Euler relation is interesting from the point of 
view of thermodynamics, because it gives a restriction on the form of state variables. 
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By the scahng laws of extensive variable and thermodynamic function, we get 

A EhiSb , A , Xb) = EbiX'^Sb , X^A , X^X^) . (4-14a) 

This denotes that Ei,{Sb, A, Xb) is the homogeneous expression of degree 1/2. Oper- 
ating the differential dx on Eq.(4-14a), we get 

^ Eb{Sb ,A,Xb) = TbSb-abA + Y,,Xb. (4-14b) 

This relation (4- 14b) is obtained from the scaling behavior (4- 14a). Furthermore by 
the well-known Euler's theorem on the homogeneous expression, the scaling behav- 
ior (4- 14a) is also obtained from the relation (4- 14b) (which is proven by the van- 
ishing differential dx[X-^Eh{X'^Sb, X'^A, X'^Xf,)] = 0). Hence Eqs.(4-14a) and (4-14b) 
are equivalent. As shown below, we find the Euler relation (4- 14b) is consistent with 
the assumption 4: 

By the Lcgendre transformation in Eq.(4-9) and the Euler relation (4- 14b), we 
get a relation, Fb = TbSb- 2abA + 2YbXb. Then substituting Eqs.(4-7), (4-lOb) 
and (4-11) into this relation, we obtain 

^1 l^ir? + ^2 iTTF + h -7. — = 2k3Xb, (4-15) 
oM on OTw 

where 

ki := -Fb (dnTb) - n {dnEb) (4- 16a) 

-r^ [{duFb) {dr^Tb) - {dr^Eb) (S^^T,) ] 

k2 := Eb {dMTb) + Tb {duEb) (4- 16b) 

-Tn^ [{drM (dMTb) " (dMEb) (dr^Tb)] 

k3 := - (dMEb) {duTb) + {duEb) {duTb) . (4-16c) 

The concrete forms of fcj (i = 1, 2, 3) are obtained from the differentials of Tb and 
Eb shown in Sec.4.1, and result in relations, 

ki = Mk3 , k2 = -Hk2. (4-17) 

Then Eq.(4-15) reduces to 

dXb dXb dXb 

This partial differential equation (PDE) is equivalent to the relation (4- 14b) which is 
also equivalent to the relation (4- 14a). Therefore, if a solution Xb of our PDE (4T8) 
exists, then the Xb satisfies the extensive scaling behavior Xb — Xb under the 
length size scaling (2-18). Indeed, the general solution of PDE (4-18) is expressed as 

Xb{M, H, r^) = rl MM/r^, Hr^) , (4-19) 

where ipb{x,y) is an arbitrary function of two arguments. Obviously this Xb is 
proportional to r^, and satisfies the extensive scaling law under the length size 
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scaling (2-18). It is also obvious that the arbitrary functions ^b{i^cfw) in Eq.(2-20) 
and ^h{.Hrw) in Eq(2-21) are consistent with ij^biM/ruj, Hruj) in Eq.(4-19), since Kc in 
Eq.(2-5) is expressed as a function of M/vyj and Hence we find that the Euler 

relation (4- 14b) is consistent with the assumption 4, which implies that the internal 
energy and also the free energy Fi, are defined well in our thermodynamics of BEH. 
The well-defined free energy guarantees the thermodynamic consistency. Now it has 
been checked that the minimal set of assumptions introduced in Sec.2 constructs the 
"consistent" thermodynamics for BEH. 



§5. Cosmological event horizon 



We examine whether the entropy-area law holds for "consistent" thermodynam- 
ics of our thermal equilibrium system for CEH. Discussion in this section goes parallel 
to Sec.4. However the integration constant in Euclidean action, which is determined 
with referring to de Sitter canonical ensemble, enables us to find a reasonable evi- 
dence of the breakdown of entropy-area law for CEH. 

5.1. Temperature and free energy of CEH 

By the assumption 3, the temperature Tc of CEH is defined by Eq.(A-4) of 
Appendix A, which relates Tc to the proper length in the imaginary time direction 
at the boundary dDc, 

-, -1 



rf3c 
Jo 



fwdr 



(5-1) 



where (3c is the imaginary time period (3-17) and is in Eq.(3-10). This Tc coincides 
with the Hawking temperature obtained by Gibbons and Hawking,^'^) and the factor 
"vZ/w is the Tolman factor which expresses the gravitational redshift affecting the 
Hawking radiation propagating from CEH to observer at r^,. Under the length size 
scaling (2-18), this temperature satisfies the extensive scaling law, Tc — Tc. 

As defined in assumption 2, the extensive state variable of system size for our 
thermal equilibrium system is the surface area of heat wall. 



A -.= 4 



7rr„ 



(5-2) 



By the assumption 3, the free energy Fc of CEH in Eq.(2-17) is defined by 
Eq.(A-5), 

Fc{Tc, A, Xc) := -Tc Iec = rcVU+ , " • (5-3) 



2 yfw 

Under the length size scaling (2-18), this free energy satisfies the scaling law of 

thermodynamic functions, F^ )^Fc. Furthermore, by the same discussion given in 
Sec.4. 1, it is also mathematically verified that the free energy Fc must be a function 
of three independent variables, which verifies the working hypothesis 1. 
Let us show some differentials for later use: 



dTc 
dM 



27rr2/, 



3/2 
cJW 



^ (3H\' -D- + ^ f 



3/^2^2 _ I 



(5-4a) 
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and 



dTc 
dH 



OF, 
dM 



AirHrjl'^ 

d[{rc- Tw) 



zmri - 1 



dF^ d[\ 



dM 



+ 



3M 3H^r^ + 1 



+ 



SH'^rj - 1 



fu 



dH 



dH 



+ 



1 



(2r^ + re - 7M - 2H\1) H 



2 2. 



(5-4b) 



(5-4c) 



(5-4d) 



where the differentials in Eqs.(2-7b) and (2-7d) are used. Then we get 



dFc _ d[{rc-ry,)^] 



dM 
dF, 

'dH 



dM 



— Trr^ 



dn 

dM ' 



c / 2 

9t Trr, 



(5-5a) 
(5-5b) 



" dH ' 

= 0, is used in the first relation. These relations 



where the definition of Vc, f{rc) 
are important to get a reasonable evidence of the breakdown of entropy-area law for 
CEH in next subsection. 

Here one might naively expect that Eqs.(5-4c) and (5-4d) would be obtained by 
replacing with rc in Eqs.(4-4c) and (4-4d), and also a relation 5m -^c = — 7i"^c ^M^c 
would be expected. However the first terms in the right-hand sides in Eqs.(5-4c), 
(5-4d) and (5-5a) appear, because of the difference of integration constant in Eu- 
clidean action as seen in Eqs.(3-8) and (3-19). The integration constant of Iec can 
not be obtained by replacing Vb with Vc in that of lEb- 

5.2. Entropy of CEH 

By the assumption 3, the entropy Sc of CEH is defined as the thermodynamic 
conjugate variable to Tc, 

dFc{Tc, A, Xc) {duFc) (dHX,) - {dnFc) (9mX,) 



dTr 



(5-6) 



where Eq.(C-3) of Appendix C is used in the second equality. Under the length size 
scaling (2-18), this entropy satisfies the extensive scaling law, Sc — >■ Sc- From this 
definition and the assumption 4 together with Eqs.(5-5), we obtain an evidence of 
the breakdown of entropy-area law for CEH by the reductive absurdity as follows: 

Assume that the entropy-area law holds for CEH, Sc = Trr^. Then Eq.(5-6) and 
Sc = irrl reduce to a PDE of Xc-, 



Jm 



dXc 
dM 



+ Jh 



dXc 
dH 



where Eq.(5-5a) is used, and 



Jm '■- 



dFc 
dH 



+ 7rr 



2dTc 
dH 



Jh :=- 



0, 



d[{rc- ru,) y/Tu 
dM 



(5-7a) 



(5-7b) 
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We find Jm ^ due to Eq.(5-5b), and Jh ^ due to dMi"c ^ and dufw = — 2/r^. 

For the case of Xc = rl^^dnbrw) given in Eq.(2-20) of assumption 4, the PDE (5-7a) 
results in a contradiction as follows: Note that, because the surface gravities and 
Kc are independent as functions of two variables (M, H) due to non-zero Wron- 
skian {duni,) {dnf^c) — {dRi^b) (dMi^c) ^ 0, the three quantities {Kb, Kc, r^) can be 
regarded as independent variables instead of {M,H,rw)- Here, the transformation 
of independent variables between two pairs {M,H,rui) and {nb, Hc,Tw) is interpreted 
as the coordinate transformation in the state space of thermal equilibrium states of 
CEH. Then we find d^^X^ = for Xc = tf^(Kfer^), and the PDE (5-7a) reduces to 



which gives Sk^^c = due to JudMi^b + JudnKb ^ 0. On the other hand, the 
form of Xc = r'^^cii^bfw) means d^f^Xc ^ 0, since ^c is not constant as explained in 
assumption 4. Hence we find the PDE (5- 7a), which is equivalent to the entropy-area 
law, contradicts Eq.(2-20) of assumption 4. 

Next, for the case of Xc = rl^dM/ru,) given in Eq.(2-21), the PDE (5 •7a) 
results in a contradiction as follows: With regarding the three quantities {Ivl , H ,Tyj) 
as independent variables, Eq.(2-21) means OhXc = and the PDE (5-7a) gives 
OmXc = due to Jm ^ 0. On the other hand, the form of Xc = r^^c{M/rw) 
means OmXc ^ 0, since 'I^c is not constant. Hence we find the PDE (5-7a), which is 
equivalent to the entropy-area law, contradicts Eq.(2-21) of assumption 4. 

The above discussions imply the breakdown of entropy- area law by the reductive 
absurdity under the minimal set of assumptions introduced in Sec. 2. Now it is con- 
cluded that we find a "reasonable" evidence of the breakdown of entropy-area law 
for CEH in SdS spacetime, where the "reasonableness" means that our discussion 
retains the "thermodynamic consistency" as shown in next subsection. 

5.3. Thermodynamic consistency of CEH 

The remaining part of Sec.5 is for the "thermodynamic consistency" of our 
thermodynamics of CEH under the minimal set of assumptions. It is also verified 
that Xc is an extensive variable and proportional to r^. The discussion given to 
BEH in Sec.4.3 is applied to CEH. 

By the assumption 3, the internal energy Ec(Sc,A,Xc) of CEH, the surface 
pressure ctc at heat wall and the intensive variable Yc conjugate to Xc are defined by 




(5-8) 



Ec 



dlnZci ^ djFc/Tc) 

diVTc) ^,^^=const. d{l/Tc) 
dFb{Tc,A,Xc) 




(5-9) 



Yc 



OA 

Eq.(4-10b) with replacing {Fb,Xb) with {Fc,Xc) 
dFc{Tc,A,Xc) 



(5-10) 



dXc 

Eq.(4-ll) with replacing {Fb,Xb) with {Fc,Xc) , 



(5-11) 
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where the relation in Eq.(5-9), Ec{Sc, A,Yc) = Fh{Tc, A,Yc) + TcSc, is regarded as 
the Legendre transformation, and Eqs.(C-2) and (C-3) of Appendix C are used in 
the second equahties in cTc and Yc. Under the length size scaling (2-18), Ec satisfies 
the scaling law of thermodynamic functions Ec ^Ec, and dc and Yc satisfy the 
intensive scaling law cTc — > cTc and Yc —>■ Yc- We find that, since the internal 
energy is a function of only extensive state variables, the state variable Xc of BEH's 
gravitational effect on CEH should be an extensive variable. (The proportionality of 
Xc to will also be shown mathematically at the end of this subsection.) 

Then by these definitions of E^ Oc and Yc together with the definition of 5c, we 
get the first law for CEH, 

dEc{Sc, A, Yc) = Tc dSc -(TcdA + Yc dXc . (5-12) 

Furthermore, the scaling behavior required in assumption 2 results in the same Euler 
relations as in Eqs.(4-14a) and (4- 14b), 

A Ec{Sc ,A,Xc) = Ec{\^Sc , \^A , X^Xc) , (5-13a) 
^ E^{Tc ,A,Xc)=TcSc-(TcA + YcXc. (5-13b) 

These two relations are mathematically equivalent by the Euler's theorem on the 
homogeneous expression. 

By the Legendre transformation in Eq.(5-9) and the Euler relation (5- 13b), we 
get a relation, Fc = Tc Sc — 2 cfc A + 2Yc Xc- Then substituting Sc-, Oc and Yc into 
this relation, we obtain a PDE of 

where l^ (i = 1, 2, 3) arc defined formally by the same definitions of ki in Eqs.(4-16a), 
(4- 16b) and (4- 16c) by replacing (F^ , T^) with {Fc, Tc). Using the differentials of 
Tc and Fc shown in Sec.5.1, we find h = M and I2 = —His which is the same 
with Eq.(417). Then our PDE (5-14) reduces to the same PDE in Eq.(4-18), and 
its general solution is 

Xc{M, H, ry,) = rl MM/r^,Hr^) , (5-15) 

where ipc{x,y) is an arbitrary function of two arguments. Obviously this Xc is 
proportional to r^, and satisfies the extensive scaling law under the length size 
scaling (2-18). It is also obvious that the arbitrary functions "^cii^bfw) in Eq.(2-20) 
and 'Fc{M/ryj) in Eq(2-21) are consistent with ipc{M /vyj^Hryj) in Eq.(5T5), since Kh 
in Eq.(2-5) is expressed as a function of M/vy, and Hvy,. Hence we find that the Euler 
relation (5- 13b) is consistent with the assumption 4, which implies that the internal 
energy Ec and also the free energy Fc arc defined well in our thermodynamics of CEH. 
The well-defined free energy guarantees the thermodynamic consistency. Now it has 
been checked that the minimal set of assumptions introduced in Sec. 2 constructs the 
"consistent" thermodynamics for CEH. 
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§6. Summary and discussions 

6.1. Summary and comments on BEH's entropy 

To research whether the thermal equihbrium is the necessary and sufficient con- 
dition to ensure the entropy-area law, we have carefuhy constructed two thermal 
equilibrium systems individuahy for BEH and CEH in SdS spacetime, and the "con- 
sistent thermodynamics" have been obtained for BEH and CEH under the minimal 
set of assumptions. The need of those assumptions was discussed in Sec. 2.1. In 
the construction of the two thermal equilibrium systems, the role of cosmological 
constant in the consistent thermodynamics has also been pointed out in the working 
hypothesis 1. In our analysis, Euclidean action method was used with referring to 
Schwarzschild and de Sitter canonical ensembles to determine the integration con- 
stants (subtraction terms). As a result, we have found a reasonable evidence for the 
breakdown of entropy-area law for CEH as shown in Sec.5.2, while the validity of the 
law for BEH could not be judged but the key issue on BEH's entropy has been clar- 
ified in Sec. 4. 2. If the breakdown of the law for BEH is verified, then it means that, 
as summarized in fourth paragraph in Sec.l, the thermal equilibrium of each horizon 
in multi-horizon spacetime is simply a necessary condition of the entropy-area law, 
and the necessary and sufficient condition of the law is the thermal equilibrium of 
the total system composed of several horizons in which the net energy flow among 
horizons disappears. 

The analysis in this paper is exact for parameter range, < -\/27 MH < 1 and 
fb < fw < fcj where the first inequality ensures that the BEH and CEH is non- 
degenerate, rf) < rc- It is obvious that our discussion is true of the near Nariai case 
(near extremal case of SdS spacetime), ~ rc ('^=^> \/27 MH ~ 1) *\ However note 
that the temperatures of horizons are equal at the exact Nariai case, Tj, = Tc at 
rb = rc- Then, Appendix D analyzes how the near Nariai case affects the entropy- 
area law. 

Turn our attention to the entropy of BEH in general case (not at Nariai limit). 
Some comments which suggest the breakdown of entropy-area law for BEH may be 
possible. Let us try to give three comments: For the first, recall the meaning of state 
variable Xf,, which expresses the thermodynamic effect on BEH due to the external 
gravitational field produced by CEH. Furthermore it is worth pointing out that the 
CEH temperature Tc depends on Kc which has dependence on (M, H), not on H 
solely. Then, in the assumption 4, it may be natural that the quantity Kc^^, not 
Hryj, is the characteristic variable of CEH's gravity and the extensive variable Xb 
is expressed by Xjj = r'^\I/b(^Hc'^w) s-s required in Eq.(2-20). If this is true, then the 
breakdown of entropy-area law for BEH is concluded as explained in Sec.4.2. 

Next, recall that, in Sec. 2. 3.1, our thermal equilibrium systems for BEH and 
CEH are compared qualitatively with the magnetized gas. By the differential of 
free energy -Fgas(Tgas, Vgas, -f^gad), the entropy 5gas and pressure Pgas of the gas are 



*' The metric of extremal SdS spacetime was found by Nariai, ' independently of the non- 
extreme SdS metric by Kottler.^^^ 
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defined by (see for example §52, 59 and 60 in Landau and Lifshitz^^^), 

(J dFgas {Tgas , Vgas ; -^ex) p _ ^-^gas (^gas i ^as ; -^ex ) /„ .. n 

'-'gas ■— ) -Tgas ■ (.D'J-j 

'^-'^gas t^'gas 

This implies that all state variables of the gas depend on the external field Hex- 
Therefore, the entropy of the gas under the influence of external magnetic field 
deviates from the entropy without external magnetic field. Hence, for our thermal 
equilibrium system for BEH, it may be naturally expected that BEH's entropy under 
the influence of external gravitational field of CEH does not satisfy the entropy-area 
law which holds for BEH in single-horizon spacetimcs. 

For the third comment, we note that, while the CEH temperature Tc is obtained 
from BEH temperature Tf, by the simple replacement of (rb , k^) with (rc , Kc), the 
CEH free energy can not be obtained from BEH free energy Fj, by such a simple 
replacement. This "asymmetry" of F^ and Fc is due to the asymmetry of integration 
constant of BEH's Euclidean action (3-8) and that of CEH's one (3-19). Then it 
is naively expected that the coefficients k of PDE (5-14) do not satisfy the same 
relation (4-17) as ki. However we find at Eq.(5-14) that the relation (4-17) holds 
for both coefficients ki and li, and the same expression of general solutions of Xf, 
and Xc are obtained as shown in Eqs.(4-19) and (5-15). This may imply that the 
same consistent structure of thermodynamics holds for BEH and CEH even though 
the forms of free energies are asymmetric. If this implication is true, then, since the 
entropy-area law breaks down for CEH, the law for BEH may also break down. 

The above three comments have no rigorous verification. Those verifications are 
left as future tasks, but the key issue is already clarified in Sec.4.2. 

6.2. Discussions 

Let us try to make two discussions. One of them is on the quantum statistics 
of underlying quantum gravity, and the other is on SdS black hole evaporation as a 

non-equilibrium process. These are independent of each other. 

For the first, we discuss about the quantum statistics: The analysis in the main 
text of this paper is based on the Euclidean action method. This is equivalent 
to assume that the basic principle of statistical mechanics of ordinary laboratory 
systems works well in calculating the partition function of the canonical ensemble 
for our two thermal equilibrium systems. Here let us emphasize that the basic 
principle of statistical mechanics of ordinary laboratory systems is the principle of 
equal a priori probabilities^'^^'^^^ *\ Hence, if the analysis in this paper and the 
comments in previous subsection are true (to imply the breakdown of entropy-area 
law for BEH) , then it suggests that the principle of equal a priori probabilities results 
in the breakdown of entropy-area law for multi-horizon spacetimes in which horizon 
temperatures are not equal and a net energy flow among horizons exits. In other 



When this principle is appUed to the micro-canonical ensemble, the Boltzmann's relation 

S = In is obtained, where S and W arc respectively the entropy and the number of states. And 
when this principle is applied to the canonical ensemble, the free energy is obtained by the relation 
in Eq.(A-5). 
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words, if the statistics of micro-states of quantum gravity obeys the principle of equal 
a priori probabilities, then the entropy-area law breaks down for the multi-horizon 
spacctimcs. 

Then what will be suggested if we adopt the other point of view? Let us dare to 
give priority to the entropy-area law, and assume that the entropy-area law holds for 
the two thermal equilibrium systems constructed in the assumption 1. Under this 
assumption, the discussion in previous paragraph implies that the principle of equal a 
priori probabilities and the Euclidean action method is not suitable to the quantum 
statistics of gravity in the multi-horizon spacetimes. In this case, the underlying 
quantum gravity should be formulated to yield the special statistic property of micro- 
states of gravity, which comes to obey the principle of equal a priori probability in 
the case of single-horizon limit. "^^^^^ 

Next turn our discussion to the second one, which is completely separated from 
the above discussion of quantum statistics. The second discussion is on SdS black 
hole evaporation process: Hereafter the heat wall introduced in the assumption 1 
is removed. Let us note the inequality > Tc due to Eq.(2-6), which means the 
existence of a net energy flow from BEH to CEH due to the exchange of Hawking 
radiation emitted by two horizons. This means that, as mentioned in Sec.l, the 
region I in SdS spacetime is in a non- equilibrium state, and the SdS spacetime evolves 
in time due to the energy flow. This time evolution is the SdS black hole evaporation 
process. Here note that Hawking temperature is usually much lower than the energies 
Eh and E,. when the horizon is not quantum but classical size.^^ Then the evolution 
of BEH and CEH during the SdS black hole evaporation can be described by the so- 
called quasi-static process, in which thermodynamic states of BEH and CEH at each 
instant of the evolution can be approximated well by thermal equilibrium states. 
(Thermodynamic state of BEH evolves on a path in the state space of thermal 
equilibrium states. Also the thermodynamic state of CEH do the same, but the path 
in state space on which CEH evolves is different from that of BEH.) This implies 
that the matter field of Hawking radiation is responsible for the non-equilibrium 
nature of SdS spacetime. Because the matter field is in a non-equilibrium state, the 
total system composed of SdS spacetime and the matter field of Hawking radiation 
is in a non-equilibrium state, even when the horizons are individually in equilibrium 
states . If we can formulate a general non-equilibrium thermodynamics for arbitrary 
matter fields which are enclosed by two thermal bodies of different temperatures, 
then a non-equilibrium SdS thermodynamics may be obtained by applying the non- 
equilibrium thermodynamics to matter fields of Hawking radiation. 

For non-self-interacting matters, a two-temperature non-equilibrium thermody- 
namics has already been constructed. Therefore, under the assumption that the 
matter field of Hawking radiation is non-self-interacting (e.g. a minimal coupling 
massless scaler field cj) satisfying D(f> = 0), the non-equilibrium evolution process 

*^ The matter field of Hawking radiation is neglected in the main context of this paper, because 
its energy scale is negligible for classical size horizons. However, when we proceed to the research 
on the non-equilibrium nature of SdS spacetime and its time evolution, it is necessary to consider 
the matter field of Hawking radiation which is responsible for the non-equilibrium nature of SdS 
evaporation process. 
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of SdS spacetime may be described by using the non-equilibrium thermodynam- 
ics of non-self-interacting matters. Indeed, the non-equilibrium thermodynamics 
of non-self-interacting matters has already been applied to the evaporation process 
of Schwarzschild black holc,^^) and has revealed the detail of evaporation process 
as a non-equilibrium process and verified the so-called generalized second law for 
the evaporation process. However the non-equilibrium thermodynamics of non- 
self-interacting matters requires to know a priori the equilibrium state variables of 
thermal bodies among which the non-equilibrium matter field is enclosed. Therefore, 
before applying the non-equilibrium thermodynamics^^ to SdS spacetime, we have 
to specify the state variables Xj, and Xc- 

Finally for self- interacting matters, under the condition that its non-equilibrium 
nature is not so strong, some non-equilibrium thermodynamics have already been 
constructed.*^^ Therefore, under the assumption that the strength of non-equilibrium 
nature is not so strong (e.g. not so large temperature difference), the non-equilibrium 
evolution process of SdS spacetime may be described by using an appropriate non- 
equilibrium thermodynamics.®-* Concerning self-interacting matters, as far as the 
author knows, no one has been applied the theories®^ even to Schwarzschild black 
hole evaporation process. 
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Appendix A 

Definitions and formulas in Euclidean action method 

Thermal field theory^^^ is a statistical mechanics for quantum matter fields in 
thermal equilibrium in flat spacetime, which is well verified experimentally for lab- 
oratory systems. It is generalized by Gibbons and Hawking^^^ to give a postulated 
formulas of temperature and free energy for thermal equilibrium systems of quantum 
gravity. This is called the Euclidean action method. Instead of the original works 
by Gibbons and Hawking, Appendix A in previous work^-* may be sufficient for 
a review of Euclidean action method. This appendix lists only the definitions and 
formulas of Euclidean action method used in the main text of this paper. 

The Wick rotation is the transformation of time coordinate t of Lorentzian space- 
time to the imaginary value —ir in the clockwise direction. 



Wick rotation -.t^—ir. 
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The Euclidean action of spacetime is defined by 

Ie[9e] := i X Iilg] with Wick rotation t ^ —ir , (A--2) 

where gE is the Euclidean metric of signature (++++) obtained from the Lorentzian 
metric g by the Wick rotation and lL[g] is the Lorentzian action of spacetime. Here 
the signature of g is { — h ++). 

In the Euclidean action method, IsigE] corresponds to the partition function 
Zci of thermal equilibrium states of the spacetime under consideration: 

In Zci := lE[gEd] , (A- 3) 

where the lower suffix means that the metric gsd the argument of Ie is the 
solution of classical Einstein equation. 

The equilibrium temperature T is defined by the proper length in imaginary 
time direction in the Euclidean space, 

-1 

(A-4) 

where /? is the period of imaginary time r determined by the regularity requirement 
of Euclidean space at the event horizon. Since the metric component gErr is a 
function of spacetime coordinates, the integral in Eq.(A-4) becomes a function of 
spatial coordinates. Therefore it is important to specify where the temperature is 
defined. Here note that the Euclidean action method is for the canonical ensemble. 
This implies the existence of a heat bath whose temperature coincides with the 
temperature of the system under consideration, since the system is in a thermal 
equilibrium with the heat bath. Therefore it is reasonable to evaluate gErr at the 
contact surface between the system and the heat bath. The contact surface is the 
boundary of the spacetime region. Hence the temperature T should be evaluated at 
the spacetime boundary. 

Then the free energy of thermal equilibrium states of the spacetime under con- 
sideration is defined by 

F := -T In Zci = -TlE[gEcl] , (A-5) 
where T is defined by Eq.(A-4) with replacing qe by gEcl- 

Appendix B 

Key points of Schwarzschild thermodynamics 

From the following key points of Schwarzschild thermodynamics, we can learn 
the minimal set of assumptions of the BEH's and CEH's thermodynamics. Indeed, 
thermodynamically consistent de Sitter canonical ensemble has already been con- 
structed with referring to the following key points in previous paper. There are 
three key points which should be shown here. The first one is the zeroth law which 
describes the existence of thermal equilibrium states: 



/ VgErrdr 
.Jo 



To what extent is the entropy- area law universal ? 



33 



Surface of Heat Bath 
(Observer is here) 



Cavity 



Fig. 3. Schematic image of thermal equihbrium of black hole with heat bath. This is described 
by the canonical ensemble. State variables of black hole are defined at the surface of heat 
bath. With those state variables, the consistent thermodynamic formulation is realized using 
the Euclidean action method.''-' 

Key Point 1 (Zeroth law of black hole) Place a black hole in a spherical cavity 
as shown in Fig. 1 and also the observer at the surface of the heat hath. Through the 
Hawking radiation by black hole and the black body radiation by heat hath, the black 
hole interacts with the heat bath. Then, by appropriately adjusting the temperature 
of heat bath, the combined system of black hole and heat hath settles down to a 
thermal equilibrium state. This equilibrium state of black hole under the contact with 
heat bath is described by the canonical ensemble. Using the canonical ensemble, the 
equilibrium state variables of black hole are defined by the quantities measured at the 
surface of heat bath where the observer is. Then the "thermodynamically consistent" 
black hole thermodynamics is obtained as summarized below. 

The second key point is the difference of black hole thermodynamics from ther- 
modynamics of ordinary laboratory systems: 

Key Point 2 (Peculiar scaling law of black hole) Extensive and intensive state 
variables of black hole show a peculiar scaling law: When a length size L (e.g. event 
horizon radius) is scaled as L ^ XL with an arbitrary scaling rate A (> 0), then the 
extensive variables X (e.g. entropy) and intensive variables Y (e.g. temperature) 
are scaled as X ^ X and Y Y , while the thermodynamic functions ^ (e.g. 

free energy) are scaled as 'P X<P. This implies that, because the system size is 
extensive, the thermodynamic system size of equilibrium system constructed in the 
key point 1 should have the areal dimension. Indeed the surface area of heat bath, 
47rr^, behaves as the consistent extensive variable of system size, where r^ is the 
radius of the surface of heat bath. 

Here recall that, in thermodynamics of ordinary laboratory systems, the in- 
tensive variables remain un-scaled under the scaling of system size, the extensive 
variables scales as the volume (for the system of three spatial dimensions), and the 
thermodynamic functions are the members of extensive variables. However, as noted 
in the key point 2, the Schwarzschild thermodynamics has the peculiar scaling law 
of state variables. Although the scaling law differs from that in thermodynamics 
of ordinary laboratory systems, the peculiar scaling law in Schwarzschild thermody- 
namics (and also de Sitter thermodynamics) retains the thermodynamic consistency 
as noted in the next key point. 

The third key point is the similarity of black hole thermodynamics with ther- 
modynamics of ordinary laboratory systems: 
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Key Point 3 (Euclidean action method and thermodynamic consistency) 

The free energy Fbh of black hole is given by the Euclidean action method, where the 
integration constant (the so-called subtraction term) of the action integral is deter- 
mined with referring to flat spacetime. The action integral is evaluated in the region 
denoted by 2M < r < ry, which is in thermal equilibrium as noted in the key point 1, 
where M is the mass parameter. Then, as for the ordinary thermodynamics, this free 
energy is expressed as a function of two independent state variables, temperature and 
system size; 

FBuiTBuATrrl) , (B-1) 

where the .rUensrve varrable Tbh := {snM^l^m^y' ^s the Hawkmg temper- 

ature measured by the observer at r^ and the factor y^l — 2M/r^ is the so-called 
Tolman factor^^^ which expresses the gravitational redshift affecting the Hawking 
radiation propagating from the black hole horizon to the observer. This Hawking 
temperature is obtained by Eq.{K-A) of Appendix A. (The temperature of heat bath 
should be adjusted to be Tbh in the key point 1.) In order to let Tbh cind 47rr^ be 
independent variables in Fbr, the mass parameter M and the heat bath radius r^ are 
regarded as two independent variables. Then the thermodynamic consistency holds 
as follows: The entropy Sbh and the "surface" pressure (Tbh are defined by 

dFBHjTBHA^rl) Abk dFBH{TBK,^7rrl) 
dT^ = ^ ' ^^«^= d^) ' ^^-^^ 

where ctbh has the dimension of force par unit area. ( See Appendix B in previous 
paper^^ for a detail explanation of the thermodynamic meaning ofaBU-) These differ- 
ential relations among the free energy, entropy and surface pressure are the same as 
obtained in thermodynamics of ordinary laboratory systems. Furthermore, as for the 
ordinary thermodynamics, the internal energy and the other thermodynamic func- 
tions are defined by the Legendre transformation of the free energy; for example the 
internal energy Ebb. i^ 

£^bh(^bh, 4Trrl) := Ebb. + Tbh Sbb. ■ (B-3) 

The enthalpy, Gibbs energy and so on are also defined by the Legendre transforma- 
tion. Then the differential relations among those thermodynamic functions and the 
other state variables also hold, for example Tbh = dEBB/dSBK- Furthermore, with 
the state variables obtained above, we can check that the first, second and third laws 
of thermodynamics hold for black holes. 

The above three key points can hold also, at least, for the other single- horizon 
spacetimes, and the consistent thermodynamics has already been established for 
those single- horizon spacetimes. ^^"^^ 

Here let us comment on the heat bath introduced in the key point 1. In York's 
consistent black hole thermodynamics,^^ the heat bath is essential to establish the 
thermodynamic consistency in the canonical ensemble as explained in the following: 
Generally in thermodynamics, as noted in the key point 3, thermod3mamic func- 
tions are defined as a function of two independent state variables. Especially the 
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free energy should be expressed as a function of the temperature and the extensive 
state variable which represents the system size. To satisfy such thermodynamic re- 
quirement, the introduction of heat bath gives us two independent variables; the 
mass parameter M and the radius of heat bath r-u,. These two independent variables 
makes it possible to define the temperature Tbh and the surface area 47rr^ as the two 
independent state variables of free energy Fbh ■ Therefore the heat bath is necessary 
to establish manifestly the thermodynamic consistency. 

Appendix C 

Useful differential formulas 







fdaj\ 




dyJ 




da2f 




\dysfj 










C.l. First case 

Let / be a function of a\, a2 and as, f = /(ai, 02,03). And consider the case 
that aj {i = 1, 2, 3) are also functions of yi, y2 and y^, ai = ciiiyj) (j = 1, 2, 3). 
Then define f{yi,y2,V's) '■= f{oii{yj)). Let us aim to express the partial deriva- 
tives df{ai,a2,cfi)/dai by the derivatives with respect to yj. Standard differential 
calculus gives, dy.f = Yli=i{9aif) i^yj^i)' which is expressed in vector form as 

dy^a2 dy^a3\ 
dy2a2 dy^a^ . (C-1) 
dy^a2 dy^asj 

Then, when det P 7^ 0, we obtain 

'da^f\ /dyj\ 

d^j\ =P-1 idyj] . (C-2) 

In section 4, by setting / = F^, (ai, 02, a^) = {Tb, A, Xb) and {yi,y2, yz) = {M, H, r^), 
we can obtain the conjugate state variables to (T^, A, Xb) as combinations of partial 
derivatives with respect to {M,H,ryj). The same is also applied in section 5. 

C.2. Second case 

Use the same definitions with previous subsection. If / has no Q!3-dependence 
(/ = f{ai,a2)) and has no ys-dependence (a^ = ai{yi,y2), i = 1, 2), then 
Eq.(C-2) reduces to 



(C-3) 



df{ai,a2) ^ {dyJ) {dy^a2) - {dyJ) {dy^aj) 
dai {dyiOLi) {dy^a2) - {dy^ai) {dy^a2) ' 

and a similar formula given by exchanging a\ and a2- 
C.3. Third case 

Use the same definitions with previous subsection. If / has no aa-dependence 
(/ = /(ai,Q;2)) and a2 has no y2- and 2/3-dependence (0:2 = Q!2(yi)) while ai 
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depends on all of yj (ai = ai (2/1,2/2,2/3) ), then Eq.(C-l) reduces to 




] (C-4) 



dy^ai 



This gives 



df{ai,a2) dy^f dyj 



dai dy^di dy^ai 

df{ai,a2) ^ jdyj) jdy^ai) - jdyj) {dy^ai) 
da2 {dyi(^2) (9^2 ai) 

^ jdyif) jdysai) - jdyj) jdy^ai) 
{dy,a2) (dy^ai) 

Furthermore, if / depends only on ai (/ = /(ai) ), then Eq.(C-5) reduces to 

dfjai) ^ dyj ^ dyj ^ dyj 
dai ^y-^^al dy^ai dy^Oii 

Appendix D 

Near Nariai case 



(C-5) 
(C-6) 



(C-7) 



This appendix analyzes the near Nariai case (near extremal case of SdS space- 
time) of our two thermal equilibrium systems of BEH and CEH. Note that, in the 
exact Nariai case, the two horizons degenerate and our two thermal equilibrium sys- 
tems of horizons disappear. Hence we consider the near Nariai case as a perturbation 
of the exact Nariai case. We introduce two independent small parameters, 5w and 
Sbc, defined by 

ritj=:3M + 5w , rc=:rb + 6bc, (D-1) 

where we require Sbc ^ M which means the near Nariai case. The parameter 5^ 
controls the position of the heat wall, and satisfies, — (3M — rb) < < rc — 3M. 

In the following, we expand the temperatures and free energies of our two thermal 
equilibrium systems by the small parameters 5^, and 5bc, in which the 0-th order 
values are of the Nariai limit dbc and 6^ 0. We measure the size of the 
exact Nariai spacetime by the mass parameter M. Then, it is useful to introduce 
the following supplemental small parameters, Sh and Sa, defined by 

H =: 5h , a=:^-Sa, (D-2) 

V27M 2 ^ ^ 



where a is given by sin a = ■\/27MH. One of three parameters (dbc, Sh , Sa) is 
independent. By the definition of a, sina = y/TTMH, and Eq.(2-3), we obtain 



2V27 MSH = 5l + 0{5i) , Sbc = 2V3MSa 



^ + lsl + o{si) 



(D-3) 
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Furthermore, from Eq.(3-10), we obtain 



36M2 



4(> 



Oi6t) + 0(616, 



(D.4) 



The requirement fw>0 means |5«;/5bc| < 1/2- 

Prom the above, we obtain the near Nariai value of temperatures (4-1) and (5-1), 



(D-5a) 



where 
Tn 



GttM 



^bc 



-1/2 



[ 1 + 0(6^) + 0{6l) + 0(6^ Si) ] . (D-5b) 



This means that, in the near Nariai case, the temperatures of our two thermal 
equilibrium systems are equal up to the leading term. Hence, at the leading term 
approximation, the total system composed of two horizons is in a thermal equilibrium 
state in the near Nariai case. Then, as discussed in the first paragraph in Sec. 6.1, 
we expect that the entropy-arca law holds at the leading term approximation in the 
near Nariai case. To see it, let us show the free energies (4-3) and (5-3) in the near 
Nariai case. 



Fb = 3M- 



3M 



Sbc 



-1/2 



Fc = Fb-3M. 



+ OiSbc) + 0{S^) + 0{Sbc (5^) (D-6a) 

(D-6b) 



Here one may think that the difference 3M = F^ — Fc results in the difference between 
entropies of BEH and CEH. But the definition of entropy, S := -dF{T, A, X)/dT, 
is important. Up to the leading term, the system size is A = 47r(3M)^ and M is 
fixed in calculating the entropy. Therefore the difference 3M does not mean the 
difference between horizon entropies. The entropy of BEH is equal to that of CEH 
at the leading term approximation in the near Nariai case. However, unfortunately, 
we can not check if the entropy-area law recovers at the leading term approximation, 
because the state variables Xi, and Xc are not specified and the partial derivative 
—dF/dA can not be calculated. At present, we simply expect that the entropy-area 
law holds at the leading term approximation in the near Nariai limit. 
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